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3 Abstract. Plasma density structures are frequently encountered in the 

4 nighttime low-latitude ionosphere by probes on the Communication/Navigation 

5 Outage Forecasting System (C/NOFS) satellite. Of particular interest to us 
e here are plasma density enhancements, which are typically observed ±15° 

? away from the magnetic equator. The low inclination of the C/NOFS satel- 
s lite offers an unprecedented opportunity to examine these structures and their 

9 associated electric fields and plasma velocities, including their field-aligned 

10 components, along an east-west trajectory. Among other observations, the 

11 data reveal a clear asymmetry in the velocity structure within and around 

12 these density enhancements. Previous observations have shown that the peak 
is change in drift velocity associated with a density enhancement occurs simul- 
i 4 taneously both perpendicular and parallel to the magnetic held, while the 

is results in this paper show that the peak change in parallel how typically oc- 
i 6 curs 25-100 km to the east of the peak perpendicular how. We discuss this 
l? and other aspects of the observations in relation to the characteristics of the 
is plasma depletions formed near the magnetic equator detected by the same 
19 probes on the C/NOFS satellite and to previous observations and theories. 
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1. Introduction 

20 Ionospheric scintillation refers to the scattering of radio waves interacting with an ir- 

21 regular ion density layer in the upper atmosphere and is an important concern for the 

22 successful operation of GPS, radar, and other communication signals. For historical rea- 

23 sons, various types of equatorial plasma density irregularities have been classified under 

24 the generic name Equatorial Spread-F (ESF) due to their common effect on radar echoes 

25 in the F-region Ionosphere [ Kelley , 2009 ]. One of the better understood irregularities 

26 associated with ESF is the plasma density depletion, also referred to as a plasma “bub- 

27 ble,” which is formed by a Rayleigh- Taylor instability acting of the bottomside of the 

28 F-layer [ e.g Fejer and Kelley 1980 ]. ESF plasma depletions appear as elongated plumes 

29 of relatively low density in radar data [ Woodman and La Hoz, 1976 ] and as curved region 

30 of reduced optical emission from space-borne imagers [ Kil et al., 2009 ]. More recently, 

31 plasma density enhancements, sometimes referred to as “blobs,” have been observed in 

32 the vicinity in plasma depletions, leading to much discussion about the potential relation 

33 between these two density structures. 

34 Plasma enhancements were first reported from ion density and temperature data from 

35 the Hinotori satellite [Oya et al., 1986 ]. They were found to occur beyond the edge of 

36 the region of occurrence for bubbles, typically about 15 degrees away from the magnetic 

37 equator in either hemisphere. Watanabe and Oya [ 1986 ] examined the statistics of oc- 

38 currence and found that enhancements favor the winter hemisphere, typically occur after 

39 midnight, and have a strong longitudinal variation pattern that varies as a function of 

40 season. Based on a statistical analysis of enhancements and depletions encountered by 
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4 1 Hinotori, they concluded that the two phenomena must have different physical drivers, 

42 with the depletions caused by ExB uplift and the enhancements caused by meridional 

43 wind instabilities. 

44 However, a number of recent observations have suggested a stronger link between the 

45 evolution of bubbles and blobs. Observations from satellites [Le et al, 2003], all-sky im- 

46 agers [ Pimenta et al., 2007; Martinis et al, 2009], radar [Yokoyama et al, 2007], and 
4 ? GPS TEC measurements [ Dashora and Pandey , 2005] have shown density enhancements 

48 forming simultaneously near depletions. Yokoyama et al. [2007] shows depletions and en- 

49 hancements forming along the same flux tube using data from ROCSAT-1 and Equatorial 
so Atmosphere Radar. Density enhancements were also found along the same held lines as 

51 ESF airglow depletions [Martinis et al., 2009]. Statistical studies with satellite data from 

52 KOMPSAT-1 and DMSP F15 have shown that the depictions and enhancements show 

53 similar seasonal and longitudinal variations (with several exceptions) [ Park et al, 2008a]. 

54 Le et al. [2003] proposes a mechanism to link the two phenomena. Polarization electric 

55 fields formed inside the depleted flux tubes are mapped down the magnetic held lines past 
se the region of depletion. These eastward polarization helds result in an upward ExB drift 

57 that transports local plasma upward to just above the depleted hux tube, resulting in a 

58 plasma density enhancement. This is supported by the similarity of the spectral properties 

59 of density fluctuations inside the enhancements and depictions [Le et al. 2003] and by the 

60 increase in the relative concentration of 0 + to H + inside the enhancement [ Yokoyama et 

6 1 al, 2007]. Recent simulation work using the SAMI3 model shows that regions of enhanced 

62 plasma density can form around spread F depletion plumes [Krall et al, 2010]. 
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es The following work discusses recent observations of plasma density enhancements from 

64 the C/NOFS satellite and their associated velocity and composition structures during a 

65 period of very low solar flux. Both plasma density depletions and enhancements were 

66 frequently encountered by the C/NOFS satellite during the first two years of operation 
e? [Haaser et ah, 2010]. While density enhancements have been observed on a number 

68 of satellites including STSAT-1, CHAMP, and KOMPSAT-1 [Park et al. , 2003, 2008b], 

69 most of these satellites were only capable of measuring density and temperature. With 

70 the exception of ROCSAT-1 observations and SAMI3 simulations, there has been little 

7 1 discussion of the ion drifts inside and around the density enhancements. The plasma drifts 

72 associated with density enhancements observed by C/NOFS will be compared the drifts 

73 associated with density depletions. 

2. Measurements 

74 The Communication/Navigation Outage Forecast System (C/NOFS) satellite is part of 

75 an Air Force mission to locate, understand, and predict equatorial ionospheric scintillation 

76 using satellite measurements, complementary ground-based observations, and physics- 

77 based models [de La Beaujardiere et al., 2004], The C/NOFS satellite was launched in 

78 April 2008 into a 13° inclination orbit with perigee near 400 km and apogee near 860 

79 km. C/NOFS consists of multiple instrument suites designed to study the ion and neutral 
so populations and their effect on the propagation of communication signals. This study will 
si focus on data from the Vector Electric Field Instrument (VEFI), provided by Goddard 

82 Space Flight Center, and the Ion Velocity Meter (IVM), provided by the University of 

83 Texas at Dallas. 
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84 VEFI consists of multiple sensors, including three electric field double probe booms, a 
as lightning detector, and a 3-axis flux-gate magnetometer [ Pfaff et al., 2010]. Of interest to 
se this study are the electric field probes and the magnetometer, which are used to provide 
s? the E xB drifts perpendicular to the magnetic field. The two perpendicular directions 
ss are defined as the meridional direction (positive outward in the plane of the magnetic 

89 meridian) and the zonal direction (positive eastward perpendicular to the meridian plane). 

90 The E xB drifts are averaged to one sample per second. 

91 IVM, part of the Coupled Ion-Neutral Dynamics Investigation (CINDI) Mission of Op- 

92 portunity on C/NOFS, consists of two sub-instruments: the retarding potential analyzer 

93 (RPA) and the ion drift meter (IDM). These two instruments provide 3D ion drifts as well 

94 as density, composition, and temperature information for the ion population [Heelis and 

95 Hanson, 1998]. The measured velocities are rotated into the field-aligned reference frame 

96 using the VEFI magnetometer measurements, providing the drift component parallel to 

97 the magnetic field. The meridional component from the IVM is also compared to that 

98 from VEFI. Unfortunately, the drift component from the RPA is not optimized due to the 

99 low ion flux during this exceptionally deep solar minimum, so the zonal component from 

100 the IVM (as well as the parallel and meridional components where the RPA contribution 

101 is significant) will be neglected for the purpose of this study. 

102 The C/NOFS data presented in this paper (Figures 1—6) are arranged into four panels 

103 per figure, each containing 5—10 minutes of data. Panel (a) shows the density ( N ) plotted 

104 on a log scale, with the total density plotted in black and the 0 + density plotted in blue. 

105 The remainder of the density consists of light ions (H + and He + ). Panel (b) shows the 

106 parallel drift (un) derived from IVM; it is defined as positive northward. Panel (c) shows 
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107 the meridional ExB drift (■ v ± m ) derived from both 1 VM (green) and VEF 1 (purple), and 
loo panel (d) shows the zonal ExB drift (v± z ) from VEFI. The important features are marked 
loo by a vertical orange line that corresponds to the peak change in density as determined by 
no a quadratic fit to the density perturbation. All data is plotted as a function of UT, with 
in the corresponding solar local time (SLT), altitude, magnetic dip angle, and geographic 
no location listed on the x-axis. 

3. Observation of Density Depletions 

ns Figures 1 and 2 are examples of plasma density depletions observed by C/NOFS. Note 
n4 that the depletions shown here are relatively weak (he., density is reduced by less than 
ns one order of magnitude) compared to other recent studies [e.g., de La Beaujardiere et a/., 
lie 2009 ] due to the difficulty in calculating accurate three-dimensional ion drifts inside the 
nr deep depletions (where density is reduced by several orders of magnitude). These “weak” 
ns depletions are similar to the upper regions of a depleted plasma shell as shown in the 
n9 simulations of Aveiro and Hysell [ 2010 ]. The two figures show that: 

120 1 . The peak change in density occurs nearly simultaneously (i.e., at the same longitude) 

121 with the peak change in velocity. This is similar to the ROCSAT observations of plasma 

122 depletions in the equatorial ionosphere [Le et al, 2003 ]. 

123 2 . A poleward change in drift is associated with most of the observed bubbles. Figure 1 

124 shows several depictions south of the dip equator. The depletions correspond to a negative 

125 (southward) change in v\\. Depletions observed north of the dip equator (Figure 2 ) have 

126 a corresponding postive (norhward) change in un. Depending on the relative strengths of 

127 the irregularity and the background drift, the direction of parallel flow may or may not 

128 change. 
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129 3. The E xB perturbation drifts within the depletions are typically upward (meridional) 

iso and westward (zonal), which correspond to a polarization electric held pointing up and 

131 east. This corresponds to the fields expected in a westward-tilted plasma bubble [Huang 

132 et al., 2010]. There are several events where the zonal drift is noisy, and the relative 

133 change associated with the bubble is unclear. 

4. Observations of Density Enhancements 

134 Figures 3-6 are examples of density enhancements observed by the C/NOFS satellite, 

135 with the associated 3D ion drifts relative to the magnetic held. A vertical line marks a 
no local density increase >50% above the background density. Plasma density enhancements 
is? are observed during solstice and equinox in both magnetic hemispheres. Figures 3-5 show 
ns density enhancements near the June solstice of 2009, with orbit tracks in the southern 

139 (winter) magnetic hemisphere (Figure 3), in the northern (summer) magnetic hemisphere 

140 (Figure 4), and near the magnetic dip equator (Figure 5). Figure 6 shows an event in the 

141 southern hemisphere near the September equinox. All four plots show similar features, as 

142 outlined below: 

143 1. The enhancements show an increase in relative 0 + concentration, displacing lighter 

144 ions. This indicates that plasma is being transported from regions with a lower apex 

145 height. This is consistent with previous observations [Yokoyama et al, 2007]. 

146 2. A poleward change in ion drift parallel to B is associated with most of the observed 

147 enhancements. Figure 3 shows several examples of enhancements in the southern magnetic 
us hemisphere; Figure 4 shows a series of events in the northern magnetic hemisphere. Note 
149 that both events are near the June solstice; consequently, the background parallel drift is 
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150 toward the southern (winter) pole. Of the nine events in these two plots, seven show a 

151 strong poleward shift in the parallel ion drift associated with the enhancement. 

152 3. The peak change in ExB drift occurs to the west of the peak change in parallel drift. 

153 The separation of the peak perturbation drifts ranges from 25—100 km along the east-west 

154 axis. The peak density is typically is found between the two peak drifts. This structure is 

155 observed in all seasons and both hemispheres; it is independent of the change in altitude 
ise of the orbit track. The apex altitudes associated with the enhancements reported here 
is? range from 483 km (Figure 5) to 1045 km (Figure 3). This structure was not reported 
us in previous in situ data. Le et al. [2003] shows that the peak density occurs at the same 
159 longitude as the peak drifts in all three directions for ROCSAT-1 data from 1999-2000. 
wo 4. These features are not limited to the higher magnetic latitudes. Figure 5 shows an 

161 event with all the characteristics of the density enhancements in the previous two figures, 

162 but it occurs near the dip equator. This may be a region of local increased density above 

163 a rising depletion as seen in simulation [Krall et al . , 2010]. 

164 5. Typically, an upward meridional drift and a westward zonal drift are associated with 

165 the density enhancements. The perturbations in the ExB drifts are consistent with those 

166 seen in the westward-tilted bubbles. A notable exception is the enhancement at 17:05:10 in 
le? Figure 6, which shows a strong eastward drift in panel (d). The perturbation drifts inside 
we this enhancement are similar to those observed inside an eastward-tilted bubble [Huang 

169 et al, 2010]. Note that both the “eastward-tilted” and “westward-tilted” enhancements 

170 observed in Figure 6 show the same longitudinal offset between the perpendicular and 

171 parallel drifts. 
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5. Discussion 

172 While the low inclination of the C/NOFS orbit restricts the longitudinal coverage of 
m the satellite at the magnetic latitudes where density enhancements are typically expected, 

174 density enhancements are frequently encountered by C/NOFS. This is consistent with 

175 the statistical studies by Park et al. [2008a], which show that the average magnetic 
we latitude where enhancements are found decreases as solar activity decreases. The following 
177 discussion is not meant to be an exhaustive study of the qualities of density enhancements 
ns during extreme solar minimum; instead it provides additional insight into the structure 
179 of density enhancements not previously discussed in the published literature. 

5.1. Comparison to previous observations 

iso The observations from C/NOFS show a number of similarities to previous observations, 

1 8 1 but are markedly different in the velocity structure across the plasma density enhance- 

1 82 rnent. As stated in the previous section, the peak change in ion drift parallel to the 

183 magnetic field occurs 25—100 km to the east of the peak change in the drift perpendic- 

184 ular to B. To better compare these events to previous in situ studies of plasma density 

185 enhancements, several of the plasma enhancements observed in data from ROCSAT-1 are 
we plotted in Figure 7. These particular events are discussed in greater detail in Le et al. 
is? [2003]. The perturbation drifts associated with the enhancement all peak simultaneously 

1 88 (i.e., at the same longitude) for the events observed by ROCSAT. 

189 It should be noted that ROCSAT has a larger orbital inclination than C/NOFS (35° 
wo as compared to 13°). Consequently, the typical orbit of C/NOFS is nearly perpendicular 

191 to the geomagnetic field lines, while the angle between a typical ROCSAT orbit and B is 

192 smaller (~50° in the case of Figure 7). Despite these differences, it is not expected that 
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193 the ROCSAT orbit track relative to the flux tube containing the density structures would 

194 mask a separation between the velocity peaks as observed in the C/NOFS data. 

195 A potential source of the difference between the C/NOFS and ROCSAT observations 
we is solar activity. The C/NOFS satellite has the distinction of being launched during the 
197 deepest solar minimum to occur since the dawn of the space age. The C/NOFS events 
we shown here were measured when F 10.7 was less than 70 sfu, while the ROCSAT event 

199 was measured when F 10 . 7 = 147 . 9 . Additionally, during 2008 the H + / 0 + transition height 

200 was significantly lower than predicted by the IRI model for an accurate input of F 10 . 7 , 

201 indicating that the ionosphere is contracted further than expected [Heelis et al, 2009 ]. 

202 Evidence of a “contracted” ionosphere is supported by the comparison between the unper- 

203 turbed regions around the enhancements in Figures 3-6 and those from ROCSAT. Despite 

204 the fact that the ROCSAT observations are 200 km higher in altitude, they are associ- 

205 ated with ion densities nearly an order of magnitude greater than those observed with 

206 C/NOFS. Figure 7 also shows a greater relative concentration of 0 + in the unperturbed 

207 background than the previous figures. 

5.2. Relation between depletions and enhancements 

208 The perturbation ExB drifts observed inside the enhancements are similar to those 

209 observed in the depletions. Both types of density irregularities typically are associated 

210 with an upward change in meridional drift and a westward change in the zonal drift. This 
on is consistent with previous studies that suggest that polarization electric fields within a 

212 density depletion are mapped down the geomagnetic field lines, forming a density en- 

213 hancement above the depleted flux tube. There is a longitudinal offset between the peak 

214 parallel and ExB drifts in the density enhancements, while the ion drifts peak in all three 
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215 directions at or near the same longitude for the depletions. The existence of a longitudi- 

216 nal offset is persistent in the observations and independent of changes in the orbit track 

217 relative to the flux tube. 

218 Simulations using the SAMI 3 model have shown that a slowly growing ESF plume can 

219 produce ExB drifts that lift density upward faster than gravity moves it downward [Krall 

220 et al, 2010 ]. This effect occurs in simulations for low solar activity (F 10.7 < 100 ) or when 

221 the growth of a plume is slowed by a mild trans-hemispheric wind. The slow growth 

222 leads to a weak super fountain [Huba et al, 2009 a] with perturbation drifts inside the 

223 irregularities on the order of 200 m/s. This is very similar to the C/NOFS observations, 

224 but the simulations do not show the longitudinal offset between the peak parallel and 

225 perpendicular drift velocities [ Krall et al, 2009 , 2010 ]. The lack of this offset may be due 

226 in part to unique conditions during the recent extreme solar minimum. 

227 The presence of a longitudinal offset may also indicate that the physics relating deplc- 

228 tions and enhancements must be better understood in three spatial dimensions. Instead of 

229 looking merely at a depleted flux tube, a three-dimensional structure is needed to explain 

230 the asymmetric velocity structure around the density enhancement. The depletion shell 

231 model described in Kil et al. [ 2009 ] would provide the necessary longitudinal asymmetry. 

232 Alternatively, zonal wind effects may also produce the necessary structure [ Huba et al, 

233 2009 b] . 

6. Summary and Future Work 

234 The ion drift velocities inside and around plasma density enhancements observed by 

235 C/NOFS are found to have a more complicated structure than previously reported. The 

236 perturbation ExB drifts associated with the enhancements are similar to those observed 
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237 within plasma depletions. However, the longitudinal separation between the peak ExB 

238 and parallel drifts is a feature previously unreported either in data or in models. It is 

239 unclear whether the longitudinal offset is unique to the recent protracted solar minimum 
2« or if it is a relatively common feature that has yet to be observed clue to the sparse coverage 

24 1 of in situ measurements of three-dimensional drifts within plasma density enhancements. 

242 This uncertainty is likely to be resolved as the C/NOFS satellite continues to take data 

243 during the transition back to solar maximum over the next few years. 

244 The plasma density enhancements in this paper were chosen to represent general obser- 

245 vations in the C/NOFS data set. A statistical study of the plasma irregularities observed 

246 by C/NOFS and the relation of three-dimensional ion drifts within these features is in 

247 progress. Further discussion of the three-dimensional nature of plasma depletion shells 

248 and associated enhancements is needed to relate these new observations to current theory. 
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Figure 1. Plasma density depictions observed south of the magnetic dipole equator. The 


panels show (a) total ion density (along with the 0 + density; the remainder is H + and He + ), (b) 


ion drift parallel to the magnetic field (positive northward), (c) ExB meridional drift (positive 


upward), and (d) ExB zonal drift (positive eastward). The three events exhibit a clear poleward 


(negative) change in the parallel drift. The apex altitude for these events is 567 km. 
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Figure 2. Plasma density depletions observed north of the magnetic dipole equator. All 


events exhibit a clear poleward (positive) change in the parallel drift. The format is the same as 


in Figure 1. The apex altitude for these events is 596 km. 
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Figure 3. Plasma density enhancements observed near June Solstice 2009 in the southern 


magnetic hemisphere. Two of the events exhibit a clear poleward (negative) change in the 


parallel drift, while the two weakest density enhancements show an equatorward (positive) drift. 


The apex altitude for these events is 1045 km. 
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Figure 4. Plasma density enhancements observed near June solstice 2009 in the northern 


magnetic hemisphere. All of the events show a poleward (positive) change in parallel flow. The 


apex altitude for these events is 636 km. 
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Figure 5. Plasma density enhancements observed near June Solstice 2009 near the magnetic 


equator. The apex altitude for these events is 483 km. 


DRAFT 


February 28, 2011, 12:30pm 


DRAFT 


KLENZING ET AL.: LOW-LATITUDE PLASMA DENSITY ENHANCEMENTS 


X - 23 


15 Sep 2009, Orbit 7675 



HH:MM:SS (UT) 17:03:00 

17:03:50 

17:04:40 

17:05:30 

17:06:20 

17:07:10 

17:08:00 

SLT (hr) 

2.0 

2.2 

2.4 

2.6 

2.8 

3.0 

3.3 

Alt (km) 

417.3 

413.4 

410.1 

407.5 

405.5 

404.2 

403.6 

Mag Dip (deg) 

-28.0 

-29.1 

-30.0 

-30.8 

-31.5 

-32.0 

-32.4 

Geo Lat (deg) 

-5.7 

-6.3 

-7.0 

-7.6 

-8.2 

-8.8 

-9.3 

Geo Lon (deg) 

132.3 

135.3 

138.3 

141.4 

144.4 

147.5 

150.6 


Figure 6. Plasma density enhancements observed near September Equinox 2009 in the southern 
magnetic hemisphere. The enhancement near 17:05:10 is associated with eastward perturbation 
drift. The apex altitude for these events is 974 km. 
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Figure 7. Plasma density enhancements observed from the ROCSAT-1 satellite. These events 


are discussed in detail in Le et al. [2003]. 
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